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Abstract 

A recently proposed strictly phenomenological static quark-antiquark poten- 
tial belonging to the generality V{r) = —Ar''°' + Kr^ + Vo is tested with heavy 
quarkonia in the context of the shifted large N-expansion method. This non- 
relativistic potential model fits the spin-averaged mass spectra of the cc, bb 
and c6 quarkonia within a few MeV and also the five experimentally known 
leptonic decay widths of the cc and bb vector states. Further, we compute 
the hyperfine splittings of the bottomonium spectrum as well as the fine and 
hyperfine splittings of the charmonium spectrum. We give predictions for not 
yet observed Be splittings. The model is then used to predict the masses of the 
remaining quarkonia and the leptonic decay widths of the two pseudoscalar 
c6 states. Our results are compared with other models to gauge the reliability 
of the predictions and point out differences. 



I. INTRODUCTION 

The charm-beauty (Be) quarkonium states provide a unique window into heavy quark dy- 
namics. The properties of the mesons are of special interest, since they are the only 
quarkonia consisting of two heavy quarks with different flavours and are also intermediate 
to the charmonium and bottomonium systems. Additionally, because they carry flavour they 
cannot annihilate into gluons so arc more stable with widths less than a hundred keV. Ex- 
cited Be meson states lying below BD (and BD* or B*D) threshold can only undergo hadronic 
(pionic) or radiative transitions to the ground pseudoscalar state which then decays weakly. 
This results in a rich spectroscopy of narrow radial and orbital excitations which are consid- 
erably more stable than their charmonium and bottomonium analogues. The discovery of 
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the Be meson by the CoUider Detector at Fermilab (CDF) Collaboration [1] in pp coUisions 
at ^/i = 1.8 TeV with an observed pseudoscalar mass Mb, (IS) = 6.40 ± 0.39 ± 0.13 GeV 
has inspired new theoretical interest in the study of the Be spectroscopy in the framework 
of heavy quarkonium theory [2,3,4,5,6,7,8,9,10,11,12,13,14,15]. 

On the other hand, for bottomonium system, the ALEPH collaboration in 2002 has 
searched for the pseudoscalar bottomonium meson, r]h, in two-photon interactions at LEP2 
with an integrated luminosity of 699 pb~^ collected at e+e~ centre-of mass energies from 181 
GeV to 209 GeV. One candidate event is found in the six-charged-particle final state and 
none in the four-charged-particle final state. The candidate /^^(IS) {r]h — > KsK~7r~^7r~TT~^) 
has reconstructed invariant mass of 9.30 ±0.02 ±0.02 GeV [16]. Theoretical estimates (from 
perturbative QCD and lattice nonrelativistic QCD) to the hyperfine mass splitting for the 
IS bottomonium state Ahf(IS) are reported (cf. [16] and references therein). 

Moreover, the inconsistent Crystal Ball measurement for the pseudoscalar charmonium 
mass, M^,(2S) = 3594 ± 5 MeV [17] has been given for more than 20 years, before a new 
measurement of 3654 ± 6 ± 8 MeV was available by Belle Collaboration in summer of 2002 
[18] for the exclusive B — )• KKsK'n'^ decays. It is close to the ?7c(2S) mass observed 
by the same group in the experiment e+e~ — >■ J/V'^c where M^,(2S) = 3622 ± 12 MeV 
was found [19]. It is giving rise to a small hyperfine splitting for the 2S charmonium state, 
AHF(2S,exp)= M^(2S) - M^,(2S) = 32 ± 14 MeV [20]. 

Very recently, several more new measurements have been appeared that support the 
Belle value: BaBar measures 3630.8 ± 3.4 ± 1.0 MeV [21] (an earlier analysis resulted in 
3632.2 ± 5.0 ± 1.8 MeV [22]), CLEO II gives 3642.7 ± 4.1 ± 4.0 MeV [23] (CLEO III prelim: 
3642.5 ± 3.6±? MeV) and a different Belle analysis yields 3630 ± 8 MeV [24]. The Ahf(2S) 
splitting calculated from a naive average of the central values of all measurements except 
Crystal Ball is 47 MeV). Further, the following mass values have been obtained: Mj/^(is) = 
3096.917 ± 0.010 ± 0.007 MeV and M^(2S) = 3686.111 ± 0.025 ± 0.009 MeV. The relative 
measurement accuracy reached 4 x 10^^ for the J/-?/'(lS), 7 x 10~® for the 'ip{2S) and is 
approximately 3 times better than that of the previous precise experiments in [25] and [26]. 
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The new result for the mass difference is M^(2S) - Mj/^(is) = 589.194 ± 0.027 ± 0.011 MeV 
[2]. 

Consequently, such observations of hyperfine splittings for the 2S charmonium and IS 
bottomonium spectra have inspired new theoretical interest in the study of the hyperfine 
splittings of the charmonium and bottomonium states as well as their spectra [5,14,27,28,29]. 
Badalian and Bakker [27] calculated the hyperfine splitting for the 2S charmonium state, 
Ahf(2S, theory) = 57 ± 8 MeV, in their recent work. Recksiegel and Sumino developed a 
new formalism [28] based on improved perturbative QCD approach to compute the fine and 
hyperfine splittings of charmonium and bottomonium [28,29]. 

As a result, one is able to test the validity of the conventional phenomenological potential 
models where the core potential follow from simple ansatzes by comparing our theoretical 
predictions of the spectrum of heavy quarkonia in terms of agreement with the experimental 
data with respect to the estimated uncertainties. Thus, these phenomenological potential 
models are not a priori connected to a fundamental QCD parameters [4,5,6,7,8]. In general, 
where the masses of bottomonium and charmonium states have been measured, the experi- 
mental uncertainties are much smaller than the theoretical uncertainties [20] . The only two 
exceptions are the very poorly measured masses of r]f,{lS) and rjc{2S). With suitable potential 
model, very good agreement with the observed spectra can be obtained for the charmonium 
and bottomonium states (e.g., [15]). These studies established the non-relativistic nature of 
the heavy quarkonium systems and, in overall, a unified shape of the inter-quark potential 
in the distance region 0.5 ^ r < 5 GeV""^. 

In this work we extend our previous analysis of the shifted large-N expansion technique 
(SLNET) developed for the Schrodinger wave equation [4, 30, 31, 32] and then apphed to 
semirelativistic and relativistic wave equations [5,32,33] to reproduce the hyperfine split- 
tings of the cc, bb and c6 spectra using a recently proposed phenomenological potential [15]. 
The calculations of the bottomonium hyperfine splittings constitute predictions of the yet 
unobserved states. The motivation of the present work is to give a detailed analysis of mass 
spectra and decay widths for cc, bb and Be systems using a recently proposed phenomenolog- 
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ical potential form [15]. We also calculate the masses of the recently found new charmonium 
?7c(2S) and the searched bottomonium ?7b(lS) mesons together with the hyperfine splittings 
of their states. 

The outhne of this paper is as following: In Section II, we first review briefly the analytic 
solution of the Schrodinger equation for unequal mass case (rriq. 7^ ^g^)- Section III is 
devoted for the spin-averaged quarkonium masses of spin triplet states. The leptonic decay 
widths are briefly reviewed in Section IV. Finally, Section V contains our conclusions. 

II. SCHRODINGER MASS SPECTRUM 

We limit our discussion to the following generality of potentials [5,15,34]: 

V{r) = -Ar-'^ + KT^ + Vo, (1) 

where A, k, a and /3 are nonnegative constants whereas Vq is taking any sign. These static 
quarkonium potentials are monotone nondecreasing, and concave functions satisfying the 
condition [4,5,30,32,33] 

V'{r) > and V"{r) < 0. (2) 

At least ten potentials of this generality, but with various values of the parameters, have 
been proposed in the literature. Cornell potential has a — P — 1, Lichtenberg potential has 
a — P — 0.75, Song-Lin potential has a — P — 0.5, and the logarithmic potential of Quigg 
and Rosner corresponds to a = — >• 0, have been recently studied (cf. e.g., [4,5,6,30,32,33] 
and references therein). The Song's potential used in [34] has a — P — 2/3. Potentials with 
a P have also been popular. Thus, Martin potential a — 0, P — 0.1 [4,5,6,30,32,33], while 
Grant, Rosner and Rynes [35] prefer a = 0.045, P = 0. Heikkila, Tornquist and Ono [36] 
tried a — 1, P — 2/3. More successful potential known in literature as Indiana potential [37] 
and the Richardson potential [38]. Recently, Motyka and Zalewski [15] have also explored 
the quality of fit in the region < a < 1.2, < /3 < 1.1 of the a, P plane. They choose the 
region with coordinates a — 1, P — 0.5. 
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In the present work, in order to get a good fit, we test the second generahty (1) with 
a = 1,P — 0.5. Therefore, the nonrelativistic phenomenological potential used by Motyka 
and Zelawiski [15] for the g^gj and qiq^ systems has the form 

^ ^ _ 0-325250 ^ Q_7Qg3g^ _ 0.78891, (3) 
r 

where V{r), ^/r and are all in units of GeV. Consequently, as r ^ 0, this potential 
has the dependence corresponding to one gluon exchange. The expected part of the 
potential which is linear in r is not seen. Probably the bottomonia are too small to reach 
sufficiently far into the asymptotic region of linear confinement. Perhaps a more flexible 
potential would exhibit the linear part, but one may be observing an effect of the expected 
screening of the interaction between the heavy quarks by the hght sea quarks [15]. The 
corresponding potential (3) is very reasonable (cf. e.g., [15]). 
The quark massese are 

rric = 1.3959 GeV, rrib = 4.8030 GeV. (4) 

For the cb quarkonium we use the reduced mass 

^ JM!}^ = 1 0816 GeV (5) 
rric + rrib 

We follow Ref.[15] in choosing the centers of gravity of the triplets for the practical reasons 
that the masses of the spin singlets for the bb quarkonium are not known or very poorly 
measured. 

For two particles system, we shall consider the N— dimensional space Schrodinger equa- 
tion for any spherically symmetric central potential V{r). If '4>{r) denotes the Schrodinger 's 
wave function, a separation of variables = Yi^jn{0,(j))u{r)/r^^~^^/'^ gives the following 
radial Schrodinger equation {h — c—1) [4,6,30,31,32] 

1 (P [^-(l-a)][^-(3-a)] ^ 



4// dr"^ 16//r 
with [I = J^^*^^ is the reduced mass for the two quarkonium composite particles. Here, En,i 
denotes the Schrodinger binding energy of meson, and k — N + 21 — a, with a representing a 



proper shift to be calculated later on and I is the angular quantum number. We follow the 
shifted 1/N or 1/k expansion method [4,5,30,31,32,33] by defining 



m=0 



(7) 



and also the energy eigenvalue expansion [4,5,30,32,33] 



oo -j^(2—m) 
En,l =^ 7^ Em-, 



m=0 



Q 



(8) 



where x — k {t/tq — 1) with tq is an arbitrary point where the Taylor's expansions is 
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being performed about and Q is a scale parameter to be set equal to k at the end of 
our calculations. Following the approach presented by Ref. [4,30,32], we give the necessary 
expressions for calculating the binding energies: 

Q 



Eo = V{ro) + 



2 ' 



(9) 



F - ^ 



Tlr + - ] U) — 



(2-«) 

SfjL 



(10) 



Eo 



Q 



16// 



(11) 



Es = %a^'\ 



(12) 



where a^^^ and a^"^^ are two useful expressions given by Imbo et al [31] and also the scale 
parameter Q is defined by the relation 



Q = 8/xroV(ro) 



(13) 



Thus, for the = 3 physical space, the Schrodinger binding energy to the third order is 
[4,30] 



En,i = V{r,) + \r,V\T,) + ^ 



;i-a)(3-a) 



k 



(2) 



(14) 
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where the shifting parameter, a, is defined by 

1/2 

(15) 

and the root, ro, is being determined via 

[8/.roV'(ro)]'/% (16) 

with = n — 1 is the radial quantum number and n is the principal quantum number. 
Once ro is found via equation (16), then the Schrodinger binding energy of the qiQj system 
in (14) becomes relatively simple and straightforward. Hence, the bound state mass of the 
qiqj system is written as 

M{qiqj)ni = ruq, + TUq. + 2E^^i. (17) 

where mq^ and m^^. arc the masses of the quark and antiquark, respectively. The expan- 
sion parameter 1/A^ or 1/k becomes smaller as I becomes larger since the parameter k is 
proportional to n which it appears in the denominator in higher-order correction. 



a = 2 - {2nr + 1) 



3 + 



V'{ro) 



1 + 2/ + {2nr + 1) 



3 + 



roV"{ro) 



III. SPIN-AVERAGED MASSES OF SPIN TRIPLET STATES 

Since the systems that we investigate in the present work are often considered as nonrel- 
ativistic system, then our treatment is based upon Schrodinger equation with a Hamiltonian 

H,^-^ + V{r) + Vss, (18) 

where Vss is the spin-dependent term. The spin dependent correction to the nonrelativistic 
Hamiltonian, which is responsible for the hyperfine splitting of the mass levels, is generally 
used in the form [7,8,39,40,41,42] 

adapted from the Breit-Fermi Hamiltonian. The number ^ substituted from the prod- 
uct of the spins corresponds to the recent assumption that the unperturbed nonrelativistic 
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Hamiltonian gives the energy of the triplet states. Since for the states with orbital angular 
momentum L > the wave function vanishes at the origin, the shift effect only the S states. 
Thus, the only first order effect of this perturbation is to shift the ^5*0 states down in energy 
by 

with the wave function at the origin is calculated by using the expectation value of the 
potential derivative via [4,5,6,30,42,43] 

In order to apply the last formula one needs the value of the wavefunction at the origin-this 
is obtained by solving the Schrodinger equation with the nonrelativistic Hamiltonian and 
the coupling constant. In this approach, the QCD strong coupling constant Q;s(4//^),on the 



renormalization point //^ is not an independent parameter. It can be connected (in the MS 
renormalization scheme) through the two-loop relation [10,15] 



2n 

^ MS 



2 



1 



In 



2fi 
MS 



(22) 



where (3q = 11 — |n/. Like most authors (cf. e.g. [4,5,6,7,8,15]), the strong coupling constant 
as{ml) is fitted to the experimental charmonium hyperfine splitting number AHF(lS,exp) = 
117 ±2 MeV [4,20] yields 

as{ml) = 0.254. (23) 

Knowing the coupling at the scale we obtain the couplings at other scales as follows. 
The number of flavours (n/) is put equal to three for 4/x^ < ml (we are not interested in the 
region 4//^ < m^), equal to four for ml > AjjF' > m^ and equal to five for 4//^ > ml (we are 
not interested in the region 4/x^ > mj). Then the value of Q;s(4/i^ = ml) from (22) is used to 
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calculate and ^\ Using the known value of A^ and the formula form (22) 

we find the value 

as{ml) = 0.200, (24) 

Prom this the value of A^~^^ is found. Note that this supports our model since a different 
choice of the Hamiltonian would in general lead to a different value of the wave function 
at the origin and to a different determination of as{m^) from the same hyperfine splitting. 
Then the estimate of as(m|) would, of course, be also different. For the hyperfine splitting 
of the c6 quarkonium we use the coupling constant 

as{4i^l,) = 0.224, (25) 

so that in each case the scale is twice the reduced mass of the quark-antiquark system. 

The calculated quarkonium masses together with hyperfine splittings are given in Tables 
1-3. The hyperfine mass splittings of cc and bb predicted by the potential model is listed 
with some other models in Table 4. Therefore, as for the hyperfine splittings in Table 
4, all of the potential model calculations try to reproduce the old experimental values, 
while the lattice calculations and perturbative QCD favor the new values. No confirmed 
experimental data to check these predictions are available yet. Let us note, however, that the 
unconfirmed experimental splitting of the 2S{cc) level —92/32 MeV- is much bigger/lower 
than expected from the potential models. In all cases, where comparison with the other 
models are significantly smaller than the splittings found by Eichten and Quigg [7] and 
similar to, but usually a little smaller than, the splittings calculated by Gupta and Johnson 
[44]. 

One can also try to compare our results with more ambitious approaches. A careful 
analysis in the framework of QCD sum rules [45] finds the hyperfine splitting of the bot- 
tomonium AHF(lS,theory)= 63^5? MeV. The central value agrees to several MeV very well 
with our expectation, but the uncertainty is too large to distinguish between the potential 
models. A lattice calculation [9] gives the hyperfine splitting Ahf (IS, theory) = 60 MeV 
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with a large uncertainty. Again the central value is close to our model, but the uncertainty 
is big enough to be consistent with all the potential models quoted here. 

IV. LEPTONIC DECAY WIDTHS 

The leading terms in the leptonic decay widths of the heavy quarkonia are proportional 
to the squares of the wave functions at the origin. Therefore, they are significant only for 
the S states. For the cc and bb quarkonium systems, we shall consider the decays of the 
n^Si (vector) states into pairs of charge conjugated charged leptons, e.g. for definiteness 
into e"^e~ pairs. For the c6 quarkonium we consider the decays of the u^Sq (pseudoscalar) 
states into TUr pairs. Since the probability of such decays contains as a factor the square of 
the lepton mass, the decays into lighter leptons are much less probable. 

The decay widths of the vector cc and bb quarkonium systems into charged lepton pairs 
are usually calculated from the QCD corrected Van Royen-Weiskopf formula [46] 

For vector mesons containing light quarks this formula leads to paradoxes (cf. [47] and 
references therein). For quarkonia, however, the main problem seems to be the QCD cor- 
rection. Thus, in order to get quantitative predictions it is necessary to include higher order 
corrections which are not known. In order to estimate the missing terms we tried two simple 
forms. Exponentialization of the first correction 

C,(a.(m^)) = exp(-i5^M), (27) 

and Padeization 

C2{as{ml)) = — (28) 

37r 

We use the average of these two estimates as our estimate of the QCD correction factor 
extended to higher orders. The difference between Ci and C2 is our crude evaluation of the 



1 - 



3n 



(26) 
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uncertainty of this estimate. The resulting leptonic widths are collected in Table 5. Further, 
we have the relation 

9 4m2 «2g2 

= SM^^^f''''^^''^' ^^^^ 

where Cav is the averaged QCD correction factor. With our choice of parameters this formula 
reduces to 

ry^-u = F{q)j^AEnF, (30) 

with F(c) = 7.07 x lO'^ and F{b) = 2.43 x 10"^ 

The formula for the leptonic widths of the pseudoscalar c6 quarkonium reads 

Tr.. = ^fk \V.,f Ms^ml (l - ^) ' , (31) 

where G is the Fermi constant, 1^6 ~ 0.04 is the element of the Cabibbo-Kobayashi- 
Masakawa matrix and the decay constant /s^ is given by the formula (cf. e.g. [48]) 

12 1^(0)1^ 2 

Be 

where C{as) is QCD correction factor. Formally this decay constant is defined in terms of 
the element of the axial weak current 

{Q\A,{Q)\BM))=^fBycbq,. (33) 

The QCD correction factor is 



C(a,) = 1 - ^ 

TT 



mb-rric, rrib 
2 m — 



(34) 



nib + n^c iT^c. 

With our parameters C{as) ~ 0.905 and since this is rather close to unity, we use it without 
trying to estimate the higher order terms. 

Substituting the numbers one finds the decay widths given in Table 5. The corresponding 
decay constants for the ground state and for the first excited S— state of the cb quarkonium 
are found to be fe, = 492 MeV and fs, = 338 MeV (cf. e.g., [6]. 
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Let us note the convenient relation 



which for our values of the parameters yields 



/gl99 ^ 

/^c = 65.2J^yA^, (36) 
where all the parameters are in suitable powers of MeV. 

V. CONCLUSIONS 

Figurel shows the behavior of the present potential form in comparison with some other 
potential forms considered in [4]. Prom this figure it is noticed that, all potential forms are 
nearly the same at large r, but the behavior seem different at very small r (r < 0.5 GeV"''^, 
about 0.1 fm), due to the variation of the one-gluon exchange term in each potential model. 
Further, the present potential model containing six free parameters: the three parameters in 
the strictly nonrelativistic phenomenological potential (3), the masses of the c and b quarks 
(4) and the strong coupling constant at the rric scale (22). This model is applicable to all 
heavy quarkonia below their strong decay thresholds. 

Consequently, we obtain the cc, bb and cb{bc) quarkonium mass spectra and also their 
leptonic decay widths in good agreement with up-to-date experimental findings. We also give 
prediction for the cb quarkonium masses and also for the leptonic widths of the pseudoscalar 
c6 quarkonium. Our model predicts similar hyperfine splitting for the IS bottomonium and 
2S charmonium as in the other potential models [7,15], lattice [49,50] and perturbation QCD 
[28,29,51]. Further, the fine splitting in charmonium is found to be M^(is) — Mj/^(is) = 597 
MeV. 

Finally, in general, the potential models reproduce the experimental values much better. 
This feature would be understandable, since the potential models contain much more input 
parameters than the lattice or perturbative QCD models. 
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FIGURES 



FIG. 1. The behavior of the different potential models versus r. 
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TABLES 

TABLE I. Schrodinger bound-state mass spectrum of cc quarkonium (in MeV). AX denotes 
the difference between the mass of particle X and the centre of gravity of the spin triplet part of 
the multiplet, where X belongs. 



State 


EQ [7] 


GJ [44] 


MZ [15] 


This work 


Expt.[20] 




3097 


3097 


3097 


3097 


3096.87±0.04 




-117 


-117 


-117 


-117 


-117 


IP(c.o.g) 


3492 


3526 


3521 


3521 


3525.3±0.2 


23^1 


3686 


3685 


3690 


3694 


3685.96±0.09 




-78 


-68 


-72 


-65.9 


-92/-32 


ID(c.o.g) 








3806 


(l3£>i state) 3769.9±25 


2P(c.o.g) 








3944 












4078 


4040±10 


2D(c.o.g) 








4150 


4159±20? 


3P(c.o.g) 








4265 












44377 


4415±6 


3L>(c.o.g) 








4430 




5^51 








4628 
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TABLE II. Schrodinger bound-state mass spectrum of bb quarkonium (in MeV). AX denotes 
the difference between the mass of particle X and the centre of gravity of the spin triplet part of 
the multiplet, where X belongs. 



State 


EQ [7] 


KR [11] 


MZ [15] 


This work 


Expt. [20] 


l^Si 


9464 




9460 


9460 


9460 


Afi-So 


-87 




-56.7 


-57.9 


(160) 


IP(c.o.g) 


9873 


9903 


9900 


9900 


9900 


2^51 


10007 




10023 


10031 


10023 


A2^So 


-44 




-28 


-23.2 




lL>(c.o.g) 


10127 


10156 


10155 


10155 




2P(c.o.g) 


10231 


10259 


10260 


10261 


10260 


33^1 


10339 




10355 


10364 


10355 


2L>(c.o.g) 




10441 


10438 


10438 




3P(c.o.g) 




10520 


10525 


10527 




43^1 








10614 




3D(c.o.g) 








10666 




53^1 








10820 
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TABLE III. Schrodinger bound-state mass spectrum of c6 (be) quarkonium in MeV. AX 
denotes the difference between the mass of particle X and the centre of gravity of the spin triplet 
part of the multiplet, where X belongs. 



State 


MZ [15] 


CK [52] 


EQ [7] 


R [53] 


G[12] 


GJ [44] 


This work 




6349 


6355 


6337 


6320 


6317 


6308 


6349 


Afi-So 


-58 


-45 


-73 


-65 


-64 


-41 


-58.2 


IP(c.o.g) 


6769 


6764 


6736 


6753 


6728 


6753 


6769 


23^1 


6921 


6917 


6899 


6900 


6902 


6886 


6926 




-33 


-27 


-43 




-35 


-33 


-30 


ID(c.o.g) 






- 








7040 


2P(c.o.g) 


7165 


7160 


7160 




7122 




7165 


33^1 














7288 


2L>(c.o.g) 














7359 


3P(c.o.g) 














7464 


43^1 














7567 


3D(c.o.g) 














7619 


53^1 














7800 
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TABLE IV. Level splittings in charmonium and bottomonium (in MeV). 



Level splitting Expt. [7]^ [IS]'' [14]^ [54]^ [49]^ [50]^ [29]'= This work 

M^{2S) - M^c{2S) 92/32 78 72 98 92 43 - 38 66 
Mt(i5) - M^^(i5) (160) 87 57 60 45 - 51 44 58 
Mt(2S) - Af^,(2S) - 44 28 30 28 - - 21 23 

'^Potential model. 
''Lattice. 

'^Perturbative QCD. 



TABLE V. Leptonic widths (in KeV). 

State EQ [7] MZ [15] This work Expt. 

1^ Slice) 8 4.5±0.5 6.72±0.49 5.3±0.4 

2^Si{cc) 3.7 1.9±0.2 2.66±0.19 2.1±0.2 

l^Siibb) 1.7 1.36±0.07 1.45±0.07 1.32±0.05 

2^51 (66) 0.8 0.59±0.03 0.52±0.02 0.52±0.03 

33^1(66) 0.6 0.40±0.02 0.35±0.02 0.48±0.08 

1^5o(c6) 4x10-8 2.8x10-8 3.58x10-^ 

2^So{cb) - 1.6x10-8 1.89x10-8 



Abb) 



{1S) = 



23 



